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A superconducting qubit coupled to an open transmission line represents an implementation of the 
spin-boson model with a broadband environment. We show that this environment can be engineered 
by introducing partial reflectors into the transmission line, allowing to shape the spectral function, 

J{u>), of the spin-boson model. The spectral function can be accessed by measuring the resonance 
fluorescence of the qubit, which provides information on both the engineered environment and the 
coupling between qubit and transmission line. The spectral function of a transmission line without 
partial reflectors is found to be Ohmic over a wide frequency range, whereas a peaked spectral density 
is found for the shaped environment. Our work lays the ground for future quantum simulations of 
other, more involved, impurity models with superconducting circuits. 

PACS numbers: 85.25.-j, 85.25.Cp, 85.25.Hv, 81.05.Xj, 78.67.Pt 


Quantum impurity models consist of a small, finite¬ 
dimensional quantum system, interacting with a larger 
system, a bath or environment. If the focus is set on 
the impurity, such as in the spin-boson model [1], these 
models constitute canonical examples of open quantum 
systems, providing insight into the validity of the Marko¬ 
vian approximation and in a variety of phase transitions 
such as the localization (i.e., suppression of the impu¬ 
rity tunneling) at large coupling strengths. However, 
such models are also extremely useful from a practical 
point of view, describing actual dynamics and coherence 
of ongoing experiments with quantum dots [2], few-level 
emitters coupled to plasmonic waveguides [3], or super¬ 
conducting qubits [4, 5]. In this context, a very relevant 
question is how to suppress or delay the qubit decoher¬ 
ence by means of external controls, [6, 7], especially as 
the coupling strengths increase. However, the focus may 
also shift from the impurity to the environment. Then, 
as in the typical case of the Kondo problem [8], we find 
that small impurities dramatically affect the transport 
properties of the bigger system. This idea has found re¬ 
newed interest in the field of quantum technologies, where 
few-level systems can be used to control photonic exci¬ 
tations, both in the optical [9, 10] and in the microwave 
regime [11, 12], or inducing photon nonlinearities [13], 
thus opening the door to all-optical integrated determin¬ 
istic quantum information processors. 

Superconducting circuits are ideal systems for the 
study of quantum impurity models and their dynami¬ 
cal properties. First, superconducting qubits are broad¬ 
band tunable devices which allow both strong and, po¬ 
tentially, ultrastrong spin-boson coupling [14-16]. Sec¬ 


ond, this coupling can be tuned in magnitude and di¬ 
rection [17-19]. Finally, state-of-the-art superconducting 
circuit technology allows the engineering of various kinds 
of photonic/bosonic environments: linear [20, 21], non¬ 
linear or even quantum [22, 23]. This flexibility permits 
a paradigm shift, where instead of dealing with the spin- 
boson model and other impurity models in complex, nu¬ 
merically difficult regimes, one can simulate these models 
by implementing them in the lab. In this way, it is possi¬ 
ble to extract both static and dynamic properties. This 
route provides access to the investigation of Kondo res¬ 
onances [24], the study of non-equilibrium phenomena 
and validation of numerical techniques [25], the testing 
of new theoretical methods for describing spin-boson dy¬ 
namics and phase transitions [26] , or even the exploration 
of quantum field theories [27, 28] . 

In this letter, we report on the experimental charac¬ 
terization of the spin-boson model using a superconduct¬ 
ing quantum two-level system (qubit) embedded into a 
broadband engineered environment, which is realized by 
a microwave transmission line containing partial reflec¬ 
tors. We report how a single qubit in a resonance fluo¬ 
rescence experiment is a powerful probe of the bosonic 
bath, providing us with both the qubit frequency Wq [25] 
and the spectral function, J(a;)jtj = i.jq- Moreover, we 
show that by introducing partial reflectors in the form 
of impedance mismatches into the transmission line, it is 
possible to shape the spectral function from the usual 
Ohmic case J{lo)(xuj to inhomogeneous distributions. 
Specifically, we study a structured environment that con¬ 
sists of two partial reflectors on each side of the qubit, 
producing a broadband, Lorentzian environment. Tun- 
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ing the qubit transition frequency over a wide frequency 
regime, we can derive the engineered spectral function, 
J{ijj) with properties intermediate between a transmis¬ 
sion line and a resonator with highly transparent mirrors 
(“bad cavity”). Finally, regarding the partial reflectors 
as point scatterers in the spirit of quantum metamateri¬ 
als [20], we also extract the spectral function of the bare 
transmission line. Our data confirms that it is Ohmic 
over a broad frequency range and allows us to extract 
the Kondo parameter. 

The two-level system in our experiment is a flux 
qubit [29, 30] built from a superconducting aluminum 
ring which is interrupted by three A1/A102,/A1 Josephson 
junctions. The effective Hamiltonian Hq = h(ujq/2) 
where is a Pauli spin operator. An external flux 
<!> allows us to tune the energy splitting of the qubit, 
hu}q = \/+ (2 Ip 6<PY. Here, the minimum energy A 
is the qubit energy gap, Ip is the qubit persistent current, 
^ — ^o/2, and is the flux quantum. The engi¬ 
neered environment is realized by a 10 mm long on-chip 
superconducting aluminum coplanar waveguide trans¬ 
mission line with a characteristic impedance Zq = 50 H, 
which is represented as Ffbath = within the 

spin-boson model. Here, a|, and Ofc are the bosonic cre¬ 
ation and annihilation operators and the oscillator 
frequencies of the bosonic degrees of freedom. The in¬ 
ductive coupling of the qubit to the engineered environ¬ 
ment is realized as a shared conductor segment of width 
500 nm and length 20 pm between the qubit loop and the 
center conductor of the coplanar transmission line. The 
corresponding interaction Hamiltonian reads as 

iJint = sin 0 ^ a\ + gl cr+ , (1) 
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Figure 1. (a) The qubit is embedded between two partial re¬ 
flectors defined by the impedance mismatch AZr = — Zq 

to the Zq = 50 line impedance. In a scattering approach, 
each segment marked by the dashed lines can be modeled 
by a transfer matrix, (b) Upper part: Calibrated transmis¬ 
sion spectroscopy data. Lower part: Zeeman shift experiment 
(data scaled by —30dB). (c) The red points: Fitted qubit 
transition frequencies. Dashed curve: Qubit hyperbola with 
A = 1.52 GHz and Ip = 630 nA. 


where are the qubit lowering/raising operators. 

The coupling rate gk of the qubit to the mode k is scaled 
by the flux-dependent quantity sin 9 = A/ujq with 9 be¬ 
ing the Bloch angle. Qubit and transmission line are 
fabricated on a silicon dioxide covered silicon wafer. The 
wafer is glued into a gold-plated copper box and mounted 
into a dilution refrigerator to a base temperature of ap¬ 
proximately 25 mK. The two interconnections between 
the on-chip transmission line and the off-chip wiring lead 
to a weak impedance mismatch at each of the two chip 
boundaries: the partial reflectors. 

The total Hamiltonian, H = Hq -\- i?bath + ^fint repre¬ 
sents an implementation of the spin-boson model with 
an engineered reservoir. In this model, the interaction 
between the flux qubit and reservoir is characterized 
by the spectral function J{uj). In the case of a bare 
transmission line without reflectors, one expects that the 
bath contribution to J(w) is proportional to w. This 
situation of a Markovian bath allows one to interpret 
^(w)|w=tjq = A(wq) as the spontaneous emission rate of 
the qubit into a quasi-one dimensional open space. In our 


experiments, we extract Ti(w) by measuring the trans¬ 
mission through the system shown in Fig. 1(a) as a func¬ 
tion of the angular frequency w, and thus can determine 
the spectral function J(w) characterizing the environ¬ 
ment. To this end, we model the transmittance through 
our system using transfer matrix theory. As detailed in 
the supplemental material [31], the total 2x2 transfer 
matrix 'T relating the incoming and outgoing modes a 
and b on one side of the system to the outgoing and in¬ 
coming modes a' and b' on the other side is given by 

T{u},U}q)= X T2uj, (uj) X Q(uj, CJq) X T2cp, (oj) X . (2) 

Here, is the transfer matrix accounting for the par¬ 
tial reflectors due to the impedance mismatches (they are 
assumed to be equal), Q describes the flux qubit and T 2 i^^ 
accounts for the propagation in the transmission line be¬ 
tween qubit and partial reflectors. We assume A^^. to 
be frequency independent and neglect potential losses in 
the partial reflectors and the transmission line. In this 
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case, the off-diagonal elements of the transfer matrix T 2 ai^ 
vanish and the diagonal ones are e^*"^ with the phase 
shift ip = 'KUjj2u}r as discussed in the supplemental mate¬ 
rial [31]. Here, ojr is the resonance frequency of the “bad 
cavity” formed in between the partial reflectors. With 
the complex transmission and reflection coefficients tz^ 
and rz^ of the partial reflectors [31], the transfer matri¬ 
ces can be written as 


Q = 
Xz, = 


1 - ^qAq 

fq/iq 

-rq/tq 

1/iq . 

ytk 


-rzjtz. 



Following Ref. 11 and Ref. 31, tq and are functions of 
the qubit transition frequency, the qubit dephasing rate 
and the qubits spontaneous emission rate Fi of the 
qubit into the open transmission line environment. As T 
is of the same form as Q, we obtain a system of equations 
for the total transmission amplitude Tt- This amplitude 
normalized to the case where the qubit is far detuned 
reads 


rp _ _ 2 (Ty i {co ^q)) 

where /\= Ti (1 — e®'^“'/“’'■)/(l-I-e®We 
note that the form of Eq. (3) suggests that the qubit 
probes a modified environment characterized by a spec¬ 
tral function J(w) = 3^(7®! (w))- In this way, Eq. (3) allows 
us to access both the spectral functions J{uj) and J{(jj) 
from a single transmission measurement as long as the 
mismatches are weak enough. 

Before analyzing the coupling of the qubit to the en¬ 
vironment in more detail, we first have to determine the 
parameters describing the flux qubit. In the upper part 
of Fig. 1(b), we show transmission spectroscopy data at 
low input power. Exemplary fits are shown in Fig. 2(a). 
In order to identify the qubit energy splitting zl, which 
lies outside our amplifier bandwidth, we use the broad 
resonance as a readout device. The results of this stan¬ 
dard Zeeman shift experiment [16] are plotted in the 
lower part of Fig. 1(b). Fitting Eq. (3) to all spectro¬ 
scopic data, we extract the parameters Z\//i = 1.52 GHz 
and Ip = 630 nA determining the flux-dependent qubit 
transition frequency uj(S<I) shown in Fig. 1(c). 

In addition to the qubit parameters, Eq. (3) also pro¬ 
vides information on the spectral function J(w) of the 
transmission line. As shown in Fig. 2(b), our data di¬ 
rectly confirms the expectation that the transmission line 
presents an Ohmic bath, J (w) / sin^ 0 = 2 tt a w, where the 
factor I/ sin^ 9 corrects for the qubit contribution to the 
frequency dependence of the spectral function. We ex¬ 
tract a Hondo parameter a ~ I % in the full frequency 
range of 3 — 7 GHz accessible by our experiments. This 
result proves the broadband nature of our experiment, 
clearly distinguishing it from standard quasi-single-mode 
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Figure 2. (a) Calibrated transmission spectroscopy data 

for a;q/27r = 3.9 GHz (top), 4.56 GHz (middle) and 5.5 GHz 
(bottom). The data points are displaced by 60dB, 35dB, 
and 10dB, respectively. Dashed lines: Fits using Eq. (3). 
(b) Spectral function J(uj) for two measurements (blue: 
Zr = 111 D; red: Zr = 253D) for the bare qubit-transmission- 
line coupling. Dashed line: Linear fit. (c) Spectral function 
J{io) for the qubit coupled to the “bad cavity”-environment 
[color code as in (b)]. 


qubit-resonator physics. The small value of a justifies our 
treatment of the transmission in the scattering approach 
in Eq. (2) and is associated with a mutual inductance 
of M~29pH between qubit and transmission line [31]. 
A mutual inductance of this magnitude corresponds to 
a single-mode coupling strength gu/^k of approximately 
15%. Such values have already enabled the experiments 
observation of ultrastrong coupling effects in systems, 
where qubits were coupled to high-quality resonators. In 
an open line, however, the quantum phase transition to 
the Hondo regime (a > 1/2) would require to increase the 
mutual inductance further to above 200 pH [31]. Next, we 
analyze J{uj) of our sample for two different, but in each 
case symmetric partial reflector configurations. These 
configurations correspond to the impedance mismatches 
Z\Zr = Zr - Zo, where Z^ = {Z, Z'} = {253 H, 111 0} and 
the associated broad resonances have center frequencies 
Wr = {a;Q,Wo} = (4.53 GHz, 4.34 GHz}. Requiring that 
the Hondo parameter of the line should not change for 
one and the same sample, we obtain the qubit gap 
A'lh = 1.2 GHz and the persistent current /{, = 644 nA. 

In the following, we discuss how an engineered reser- 
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Figure 3. (a) Uncalibrated transmission power with the qubit 
far detuned (diamonds). Dashed line: Lorentzian fit of the 
resonance, (b) Lamb shift of the “bad cavity”-resonance. 
Dashed line: Fit using Eq. (4). 


voir, as we apply it in our experiments, can be used to 
alter the frequency dependent coupling of the qubit to 
the transmission line. J(w) is comprised of a resonator- 
type contribution, a Lorentzian function characterized 
by the loss rate k, and an Ohmic line contribution 
= 2'Kauj characterized by the Kondo parameter a. 
Our model involves both of these contributions. Con¬ 
sequentially, the radiative correction factor applied to 
the spectral function of the Ohmic line tends towards 
the Purcell factor for k—^ 0 [31]. In contrast to this ex¬ 
treme limit, however, we conduct our experiments in a 
regime where both contributions matter. In this way, 
we are able to extract the spectral function of the un¬ 
derlying Ohmic line J(w) in a broad frequency range. 
In Fig. 2(c), we show that we can, by varying the re¬ 
flectance of the partial reflectors, effectively shape the 
total spectral function J{uj). On the center frequency 
of the Lorentzian function, J(a;o)/2 tt sin^ 0 ~ 8.25 GHz 
and J(a;Q)/2 TT sin^ 6* ~ 2.33 GHz. Furthermore, from 
the fits of Eq. (3), we obtain Ktra/27r~ 185MHz and 
756 MHz as the decay rates for the resonances 
connected to the center frequencies Wq and Wq, respec¬ 
tively. 

Next, since our experiments happen in a regime inter¬ 
mediate between resonator and transmission line physics, 
we can also analyze our data from a resonator-qubit point 
of view. We are, however, limited by the quality of the 
Lorentzian fits and therefore only show results for the 
measurement with Zi. = 25317. For this reflector config¬ 
uration, Fig. 3(a) shows a Lorentzian fit to the uncal¬ 
ibrated transmission spectrum. The decay rate of the 
resonator, Kj- = 169 MHz, in the case that the qubit is 
far detuned agrees well with Ktra found from the trans¬ 
mission line analysis. Furthermore, the center frequency 


of the Lorentzian, Wo,L) varies as a function of the flux 
threading the qubit. As shown in Fig. 3(b), this Lamb 
shift follows the expected behavior [32], 

wo.L=wo+5^ (sin6»)^ (--^-- ) . (4) 

yWq—Wr Wq-t-Wr/ 

A numerical fit to the data yields a coupling rate 
5 sin 0 = 216 MHz. To compare the transmission line 
coupling with g, we calculate (c.f. supplemental mate¬ 
rial [31]) a coupling gtra = 233 MHz the qubit would have 
to a single mode resonator with resonance frequency Wr 
made from our transmission line. Remarkably, we find 
a very good agreement between the value of g extracted 
from a resonator type experiment with the value of gtra 
extracted from the qubit fluorescence. 

Finally, we discuss the increase in the spontaneous 
emission rate of the qubit for our engineered environ¬ 
ments. In this context, we remark that the expres¬ 
sion J{uj) = ^{ri{uj)) only holds for Markovian envi¬ 
ronments where excitations emitted to the environment 
are not allowed to excite the qubit again. Obviously, 
the latter assertion breaks down with the onset of the 
strong coupling regime. Using this, we can by setting 
«^tra(wr) = 5tra(wr) define a Critical decay rate /«crit(<^r)- 
With /«crit(wo) = 233 MHz and Kcrit(wQ) = 175 MHz, we 
find that in the first case the Markovian approximation 
of the bath is violated (wtra ^ Kcrit(<j-’o)). This situation is 
different in the second case, where ^ Kcrit(wo), and 
we can interpret A/A = J(loq)/ J(ujq) 8 in terms of an 
increased spontaneous emission. 

In summary, we experimentally implement the spin- 
boson model with an engineered reservoir using super¬ 
conducting circuits. To this end, we study a supercon¬ 
ducting flux qubit in a continuous environment modified 
by two discrete scatterers. The formation of a broad reso¬ 
nance effectively allows us to control the density of states 
in the environment of the qubit. In this way, we can in¬ 
crease the spontaneous emission by a factor of 8. Our 
model of the combined qubit-resonance system addition¬ 
ally supplies information about the spontaneous emission 
rate of the qubit in a bare transmission line. We confirm 
that the transmission line represents an Ohmic environ¬ 
ment and our coupling corresponds to a Kondo parameter 
a = l%. Additionally, we crosscheck our model by com¬ 
paring the resonator Lamb shift for the qubit being far 
detuned and find a very good agreement to the transmis¬ 
sion line coupling. Our work constitutes an important 
step towards the quantum simulation of impurity mod¬ 
els with superconducting circuits. In this context, the 
modeling of the effect of broad resonances on the envi¬ 
ronment paves the way for protocols where qubits can 
be coupled to and decoupled from engineered broadband 
reservoirs dynamically in the spirit of band engineering 
in photonic crystals. Finally, we show that increasing the 
mutual inductance to the technologically challenging but 
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in principle feasible value of 200 pH is required to reach 
the Kondo regime with flux qubits. 
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I. EXPERIMENTAL TECHNIQUES 

All measurements are conducted in the dilution refrigerator at a base temperature of approximately 25 mK. A 
schematic diagram of the experimental setup with the sample under study and the different circuit elements placed 
at different temperature levels is shown in Fig. 1. A large magnetic shielding is obtained by three mu-metal shields 
at room temperature, a cryoperm shield around the 4K stage and a superconducting lead shield thermally anchored 
to the still (0.7 K). The flux bias for the qubit is created by the field of an external coil which can be operated 
in the persistent current mode and hence applies minimal flux noise to the sample. We measure the microwave 
transmission through the sample using a vector network analyzer. The input line is heavily attenuated at various 
temperature stages to minimize thermal noise at the sample. At the output port a cryogenic HEMT amplifier and 
a room temperature RF amplifier are used to amplify the signal while two isolators protect the sample from the 
amplifier noise. For two-tone experiments, the microwave transmission through the sample is measured by the vector 
network analyzer while a second microwave tone is supplied to the input line by an additional microwave source. 




Figure 1. Experimental setup for measuring the microwave transmission through the sample chip in a dilution refrigerator with 
base temperature of 25 mK using a network analyzer. The red boxes display attenuators. 


II. TRANSMITTANCE OF AN ATOM EMBEDDED INTO A ID TRANSMISSION LINE WITH TWO 

PARTIAL REFLECTORS 

In this section we present a simple model that can be used to analyze the transmittance through a one-dimensional 
transmission line containing an atom which is surrounded by two partial reflectors. The latter originate from small 
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impedance mismatches between the on-chip superconducting coplanar waveguide transmission line and the off-chip 
wiring. Our model is based on the language of transfer matrices and fully reminiscent of scattering theory in one¬ 
dimensional systems. We first address the issue of a pseudo-resonator that is created by the two impedance mismatches 
and compute the corresponding bound states. The resonant structure related to this resonator is referred to as weak 
resonance in the main text. In the next step we couple the resonator to a qubit which is placed in the center of this 
resonator and acts as a scatterer. We then calculate the new transmission and reflection properties of the hybridized 
system. In particular, we show that the qubit is affected by the change in the density of states (DOS) induced by the 
weak resonator. Compared to the ideal transmission line (open space scenario) the effect of the weak resonator can 
be accounted for by a radiative correction. 


A. Transmission and Reflection without Qubit 

We consider a superconducting transmission line with impedance Zg = 50 Q containing two partial reflectors induced 
by small impedance mismatches AZ-^ = Z,- — Zq at the transition to the off-chip microwave wiring. These partial 
reflectors result in a finite scattering of the propagating microwaves and can be modeled by scattering or transfer 
matrices. We note that a scattering matrix usually connects the incoming modes a and a' on both sides of the 
scatterer to the outgoing modes b and b', whereas in the transfer matrix formalism the transfer matrix relates the 
incoming and outgoing modes a and b on one side of the scatterer to the outgoing and incoming modes a' and 6' 
on the other side. In the following we will stick to the transfer matrix formalism. The transfer matrix captures the 
dispersive properties of the partial reflectors and can be written explicitly in terms of frequency-dependent reflection 
and transmission coefficients. 

The general form of a transfer matrix for a partial reflector can be written as^’^ 





rz, 
\ tz^ 



1 


tz, / 


G SU(1,1) , 


( 1 ) 


where tz^ and rz^ = AZ,-/(Zq + Z^) are the complex transmission and reflection coefficient, respectively. If the scat¬ 
tering process is lossless, \tz^\ + |cz,.| =1. The transfer matrix relates the incoming and ontgoing wave amplitudes 
a{b') and b{a') on both sides. 


b' 



( 2 ) 


where from now on we strictly keep a right moving direction for the transfer matrices. 

Note that G SU(1,1) is written in a fully symmetric form^ with time reversal symmetry and identical trans¬ 
mission and reflection coefficients for waves incident from either side. The unitarity of the corresponding scattering 
matrix^ implies conditions, i.e. tz^ + rz^ = 0, for the orthogonality of the respective eigenmodes. In particular, 

the symmetric reflectivity coefficient is always 7r-shifted from transmittance, i.e. rz^/tz^ = i x K. Taking into account 
= 1, the reflectivity 


rz. 




(3) 


can always be specified in terms of transmittance. 

Our model considers a setup of equal partial reflectors connected by a 50 Q transmission line forming a pseudo¬ 
resonator. The transfer matrix C describing this pseudo-resonator is given by 


n = Xz^-T^^-Xz, , (4) 

where we use 

iks Q 

0 e-''=® 



( 5 ) 
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in order to account for the free and lossless propagation of plane waves in between the two partial reflectors separated 
by a distance s. 

The distance s between the partial reflectors can be measured in units of the resonant mode wavelength Ar of the 
weak resonator. For the fundamental mode we have s = Ar/ 2 . Then, with k = 2 tt/\ and a linear dispersion uj = ck we 
can express the optical path fc • s as 


Aj. UJ 

k ■ S = TT— = TT-. 

A UJr 


( 6 ) 


For the transfer matrix TZ we then obtain the well known result of the transmission through a double tunneling 
barrier: 


7^ = 


/ ^ _R^\ 

T,* Tr* 
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V Tr J 
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where 


Tr = 


^iks J.2 
^ ^Zr 


1 — 


, „ tz. rz,-r*z^ 

and Kr = —— 


1 _o2iA;s 

^ ^ Zr 


We immediately see that reflection and transmission are 7 r-phase-shifted, i.e., 

i?r _ 2 - Im (e*''® rz^) 


Tr 


1*2 


( 8 ) 


(9) 


This implies Rt*T,- + R^- T^* = 0 and allows us to introduce normalizable bound states for the pseudo-resonator. Con¬ 
sidering waves incident from opposite directions, we seek two scattered waves that are orthogonal, i.e.^ 


/ OO 

da; (l>T*{x) 01 (a;) = 0 . 

-OO 


( 10 ) 


If we place the origin a; = 0 at the center of the resonator, we find for the left incoming wave with uniform scaling 


L: 


01 (a;) = 




^ik{x+s/2) ^-ik{x+s/2) 


Il_ Qik{x-s/2) ^-ik{x-s/2)^ 

tz^ ' tz^ 


and similar for the right incoming wave: 


Xj. Q^k(x-s/2) 
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( 12 ) 


ei'=(^-®/2) + Q-ik{x-s/2)^ ^ ^ 


The latter can be formulated in terms of the symmetric and anti-symmetric superpositions, i.e. 

0±(a:) = [01 (a;) ± 0r(a:)] . 


(13) 


We are only interested in the intra-resonator solutions ' 0 ±( 2 :) = 0±(a:) with —s/2 < x < s/2. From Eqs. (11, 12) 
we hence obtain 


0_(a:) =i^—fe-*'=®/^-e''=®/^r 2 ,)sin(fca:) , 0+(a;) = ^ -f cos(fca;) , (14) 


1 T, 
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which yields 


\ip-{x)\ = 


1 + rz^ 


and |V.+ (a)P = 


cos^(a;) 


L 


1 — e' 


ifcs 




(15) 


B. Transmission and Reflection with Qubit 


The most general form for a transfer matrix with time reversal symmetry and symmetric reflectance and trans¬ 
mittance tq can be written as^ 


Q = 



V tq 




(16) 


We note that the general transfer matrix Q may have lower symmetry than the transfer matrix in Eq. (1) if it 
does not conserve probabilities, i.e.,|tq| -I- |rq| ^ 1. This allows us to consider possible dephasing effects. This is 
particularly useful when the transfer matrix Q is applied to to model the scattering event with a qubit. 

In order to obtain the transfer matrix T for coupled qubit-resonator-system, we extend Eq. (4) and place the qubit 
in the centre of the resonator. The total transfer matrix then reads as 


T = Xz^ ■T2 ui^ - Q ■T2ui^ ■ Xz^ ■ 


(17) 


Without specifying transmittance tq and reflectance rq for the qubit, we obtain 


with 


and 


T = 
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Tt 




Tt = tq 


1 - 2rqrz^ 


e*'=" tl 



02ifes 


tz^ (l + kzj") 

t*^^ 1 — 2rqrz^ — r'^ (tq^ — Xq^) e^‘^® 


(18) 


(19) 


( 20 ) 


The remaining task is to derive expression for the transmission and reflection coefficient of the qubit. We start 
from Ref. 5 and assume low power and a unit coupling efficiency between qubit and the line field. Then, for a qubit 
that emits with the same phase in both directions (tq — Tq = 1), we have for any arbitrary relaxation rate Fi and 
dephasing rate F^p 




2F^-iA 


( 21 ) 


with A = oj — ojq being the detuning from the qubit transition frequency Wq. 

Note that the transmission coefhcient given by Eq. (21) is applicable to describe the resonance fluorescence of the 
qubit in open space. In this case we can assign a physical meaning to Fi and Fp,: Fi is the spontaneous relaxation 
rate and the pure dephasing rate®. 

In order to show that this meaning can be restored also for a transmission line containing two partial reflectors, we 
calibrate the transmittance in Eq. (19) for the hybridized system with respect to the transmittance of the qubit-free 
case in Eq. (8). For the calibrated transmittance Tcai, we obtain 


1 riz j, ^ 2{F^-i{uj-u}q)) 

2Fp-iA FiZ + 2Fp-2i{uj-ujq) 


( 22 ) 
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with 


l + e*'=" 

1 — e*''® rz^ 


(23) 


Evidently, the transmission coefficient of the coupled qubit-resonator-system is obtained from the one for a qubit in 
open space by multiplying with a radiative correction factor. 

With Eq. (9) and Eq. (15) we can rewrite the latter in terms of its real and imaginary part, i.e. 


V TJ |^^( 0 )| 


(24) 



Figure 2. Amplification factor of the spontaneous emission rate due to Purcell effect (dashed line) or our model (diamonds). 

We stress that the qubit couples only to the symmetric cavity mode (a;) since the anti-symmetric mode vanishes 
completely at the position of the qubit, i.e. '0_(O)=O. The real part of 2 is consistent with the calibrated DOS at 
the position of the qubit inside the cavity, i.e. ip+{^) = with t —)■ 1. With i?c/T'c = i x K, the imaginary part of 

z marks a Lamb-shift for the resonance of the calibrated transmittance Tcai- 

In the remainder of this section, we compare our model to the Purcell effect where the density of states of the 
bath is assumed to be resonator-like. The Purcell factor measures the change in the spontaneous emission due to the 
altered density of states. The one-dimensional representation of this factor reads^ 

Xpurcell <2 i ( /^^2 ’ 

where k is the decay rate of the resonator. In Fig. 2, we compare our model to Eq. (25). As expected, for small n 
values our model tends towards the Purcell factor. 


III. SPIN-BOSON MODEL OF A QUBIT IN AN OPEN TRANSMISSION LINE 

We start with the interaction Hamiltonian of a flux qubit coupled to a transmission line® 

iLint = (skOk +5kak) ■ (26) 

k 

Here, ax is a Pauli spin operator and is the coupling to the field mode (ok -I- aj^). On the other hand, for an 
inductive coupling we can write 


Hint — Af /line i 


(27) 
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where M is the mutual inductance and /une is the vacuum current of the line. The vacuum current reads as 

1 


.^line — ^ ^ ^ 




(28) 


Here, Iq is the self-inductance per unit length, Z is the line impedance and L is the mode volume. Using Eq. (26-28), 
we obtain 



with 

G = iMIp^ 



(29) 


(30) 


In the next step, we want to model a situation where the qubit only couples to one specihc mode Wk equivalent to 
resonator coupling. In closed boundary conditions 


5k = G 



sin(A:a;) 


(31) 


with the normalization L = A/4. 

In the spin-boson model, the spectral function J{u}) describing a quantum two-level system coupled to the environ¬ 
ment consisting of an infinite number of bosonic degrees of freedom is defined as 


J{uj) = 2 tt ^5k5k - ‘^k) ■ 

k 


(32) 


Using Eq. (31) and moving to a continuum description, we find 

GG* 

J(a;) = —D(u;)w , (33) 

V 

where v is the phase velocity of the wave and D(uj) is the density of states in the line. In order to compare the 
transmission line coupling with the resonator coupling, we assume that the transmission line couples only to a single 
mode. Using Dioj) = 2, we find 


l5tra| = . (34) 

Sometimes it is convenient to express the Ohmic spectral function in terms of the dimensionless Kondo parameter* 
J{uj) = 2'Kauj. From Eq. (33) and Eq. (30) we can derive 

GG* 

a=- - D{lu)(xM^ . (35) 

2 TT U 
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